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Abstract 

A global event shape analysis of the multihadronic final states observed in neutral 
current deep inelastic scattering events with a large rapidity gap with respect to the proton 
direction is presented. The analysis is performed in the range 

5 < Q2 < GeV^ and 

160 <W< 250 GeV, where is the virtuality of the photon and W is the virtual-photon 
proton centre of mass energy. Particular emphasis is placed on the dependence of the shape 
variables, measured in the 7*— pomeron rest frame, on the mass of the hadronic final state, 
Mx- With increasing Mx the multihadronic final state becomes more collimated and 
planar. The experimental results are compared with several models which attempt to 
describe diffractive events. The broadening effects exhibited by the data require in these 
models a significant gluon component of the pomeron. 
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1 Introduction 



At HERA a class of neutral current (NC) deep inelastic scattering (DIS) events has been ob- 
served which is characterized by a large rapidity gap (LRG) between the proton beam direction 
and a multi-particle final state [||, H. The dependence of the LRG event fraction points to 
a leading twist production mechanism. The rate of LRG events is most commonly accounted 
for by introducing "diffractive" processes, which proceed via the t-channel exchange of a colour 
singlet object with quantum numbers of the vacuum, called the pomeron (iP). However, the 
nature of the pomeron is at present far from clear. Ingelman and Schlein assumed that the 
pomeron emitted from the proton behaves like a hadron and suggested that it could have a 
partonic substructure which could be probed by a hard scattering process. The UA8 experi- 
ment at CERN later observed events containing two high-py jets in interactions which were 
tagged with leading protons (or antiprotons). This observation was explained in terms of a 
partonic structure in the pomeron [^. 

Jet production in LRG events at HERA in the photoproduction regime and the pattern of 
scaling violations in diffractive DIS led to the conclusion that in order to interpret the data 
within the Ingelman- Schlein approach, the pomeron must have a high gluon content. 
Various models for diffractive scattering can be used to describe the hadronic final state in 
LRG events [0. The theoretical predictions range from a calculation of gluon pair exchange 
P, H, to a phenomenological description of pomeron quark and gluon densities [|10], |ll|] and to 
a boson-gluon fusion scheme, where the resulting quark pair evolves into a colour singlet state 



121 , |T^]. Recently, calculations for two- and three-jet production in DIS rapidity gap events 
have been presented, in which a scalar pomeron is considered to have a pointlike coupling to 
quark and gluon pairs [1^ . 

In order to obtain a deeper understanding of multihadronic final states with a large rapidity 
gap, we study shape variables in the centre of mass frame of the observed hadronic system (of 
mass Mx)i interpreting this as the ^*—lP cm. system, and investigate their dependence upon 
Mx- This is in analogy with the studies of shape variables in e^e~ annihilation as a function 
of \fs [0 and to those which led to the interpretation of three-jet events as a consequence 
of gluon bremsstrahlung |T^, |T^. It has also been applied to the dependence of fixed 



target antineutrino data as a function of W 



2 Experimental setup 

The data presented here, taken with the ZEUS detector at HERA in 1994, correspond to a 
luminosity of 2.57±0.04 pb~^. HERA operated with 153 colliding bunches of 820 GeV protons 
and 27.5 GeV positrons. Additional unpaired positron and proton bunches circulated, which 
were used to determine beam related backgrounds. 

A description of the ZEUS detector can be found in P2| , PT[| . The components used in 
this analysis are briefiy discussed here. The uranium-scintillator calorimeter (GAL) p3[ covers 
99.7% of the total solid angle. It consists of the forward calorimeter (FCAL) covering the 
range in pseudorapidityQ 4.3 < r] < 1.1, the barrel calorimeter (BGAL) covering 1.1 < t] < 

^The ZEUS coordinate system is defined as right-handed with the Z axis pointing in the proton beam 
direction, and the X axis horizontal, pointing towards the centre of HERA. The pseudorapidity is defined as 
rj = — ln(tan |), where 6 is the polar angle with respect to the proton direction. 
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—0.75 and the rear calorimeter (RCAL) covering —0.75 < rj < —3.8. Each calorimeter part is 
segmented into electromagnetic (EMC) and hadronic (HAC) sections. Each section is further 
subdivided into cells of typically 5 x 20 cm^ (10 x 20 cm^ in the RCAL) for the EMC and 
20 X 20 cm^ for the HAC sections. Under test beam conditions the calorimeter has an energy 
resolution of a/E = 18%/-^ E{GeV) for electrons and a/E = 35%/-y^(GeVy for hadrons. 
The timing resolution of a calorimeter cell is less than 1 ns for energy deposits greater than 
4.5 GeV. In order to minimize the effects of noise due to the uranium radioactivity on the 
energy measurements all EMC(HAC) cells with an energy deposit of less than 60(110) MeV 
are discarded from the analysis. 

The tracking system consists of a vertex detector (VXD) a central tracking chamber 
(CTD) p5|, and a rear tracking detector (RTD) p6|, ^ enclosed in a 1.43 T solenoidal magnetic 
field. The interaction vertex is measured with a typical resolution along and transverse to the 
beam direction of 0.4 cm and 0.1 cm respectively. 

The position of positrons scattered at small angles with respect to the positron beam direc- 
tion is measured using the small angle rear tracking detector (SRTD) which is attached to the 
front face of the RCAL. The SRTD |^ consists of two planes of scintillator strips, each one 
1 cm wide and 0.5 cm thick, arranged in orthogonal directions and read out via optical fibers 
and photo-multiplier tubes. 

The luminosity is measured via the Bethe-Heitler process, ep e'yp, using a lead-scintillator 
calorimeter (LUMI) which accepts photons at angles < 0.5 mrad with respect to the 
positron beam direction and is located at Z = —107 m. 



3 Reconstruction and kinematic variables 

The kinematic variables used to describe deep inelastic ep scattering 

e{k) + p (P) e {k') + X 

are the following: 

the negative squared four-momentum transfer carried by the virtual photon 

Q' = -q' = -{k - k')\ 

the Bjorken x variable 

X 



2P-q 

the fractional energy transfer to the hadronic final state 

P -q 

y = T-k^ 

and the square of the centre-of-mass energy of the virtual-photon proton system (7*p) 

= {q + Pf = ^^^^^—^ + Ml 
X ^ 

2 



where Mp is the proton mass. 

These variables, only two of which are independent, can be determined either from the 
scattered positron or from the hadronic system. They can also be determined from a mixed 
set of variables, the so-called double angle method . In this analysis, the scattered positron 
method is used because it provides a better resolution in the kinematic range under study. 
The variables and calculated from the scattered positron variables, are given by the 
expressions 

Ql = 2E,E'^{l + cose'^), 

and 

2/e = 1-^(1- CoO, 

where is the incident positron energy and E'^ and 9'^ denote the energy and angle, with 
respect to the proton direction, of the scattered positron. 

To identify the neutral current DIS events the quantity 5 is also used: 

5 = ^{Ei-pzi), 
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where Ei denotes the energy assigned to the signal in calorimeter cell i and pzi = Ei cos 6., 
9i being the polar angle of the cell. The sum is running over all calorimeter cells. For fully 
contained DIS events, and neglecting detector smearing effects, 5 = 2- E^. 

In order to describe the kinematics of diffractive scattering 

e {k)+p {P)^e {k')+p {P')+X 

the following two additional variables are introduced: the squared four-momentum transfer at 
the proton vertex, 

t = {P- P')\ 

whose absolute magnitude is small compared to + M|- in the diffractive processes selected 
in this analysis, and Xjp, the fraction of the proton momentum carried by the pomeron, defined 
as 

^ (P -P')-q ^ Ml + Q^-t ^ Ml + Q'' 
~ P-q ~ W^ + Q^-M^~ W^ + Q''' 

where Mx is the mass of the hadronic system. 
In addition, the variable (3 is defined as 

X Q2 



P 



2{P-P')-q xr Ml + Q'^-t Mj. + Q^' 
which can be interpreted as the fraction of the pomeron momentum carried by the struck quark. 
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4 Trigger and data selection criteria 



Events were filtered online by a three level trigger system ^T[|. At the first level DIS events 



are selected by requiring a minimum energy deposition in the EMC section of the CAL. The 
threshold varies between 3.4 GeV (in the RCAL) and 4.8 GeV (in the BCAL). At the second 
level (SLT), non-ep background is further reduced by using the measured times of energy 
deposits and the summed energies from the calorimeter. The events are accepted if 
5sLT = Ei{l — cos9i) > (24 — 2Ej) GeV, where Ei and 9i are the energies and polar angles 

i 

of calorimeter cells, and is the energy deposit measured in the LUMI photon calorimeter. 
This requirement on Sslt does not reject events with initial state QED radiation. 

The full event information is available at the third level trigger (TLT). Tighter timing cuts as 
well as algorithms to remove beam halo and cosmic muons are applied. The quantity 6tlt was 
determined in the same way as 6slt- The events were required to have 6tlt > (25 — 2Ey) GeV. 
Finally, events were accepted if a scattered positron candidate of energy greater than 4 GeV 
was found. In total 900,853 NC DIS candidates satisfied the above trigger conditions. 

DIS events were selected off-line by looking for scattered positron candidates recognized using 
the pattern of energy deposition in the calorimeter. The efficiency of the positron identification 
algorithm, as well as the procedures to incorporate information from the SRTD in order to 
improve the positron energy and scattering angle measurements are thoroughly discussed in 
our recent publication on the measurement of the proton structure function F2 [pT| . 

The off-line cuts used to select multihadronic DIS events are the following: 



> 8 GeV. 
5 GeV^ <Ql< 185 GeVl 
l/e < 0.95. This cut removes fake positrons found in the FCAL. 

35 GeV < 5 < 65 GeV. This cut removes events with large initial state QED radiation 
and further reduces the background from photoproduction. 

The events are required to have a reconstructed vertex. This cut suppresses beam-gas 
background. 

Events with a scattered positron impact point in the RCAL inside a box of 26 cmx26 cm 
around the beam pipe are rejected. This guarantees full containment of the shower in the 
calorimeter. 



A total of 294,527 events were selected this way. 

The mass of the final state hadronic system, Mx, can be reconstructed from the energy 
deposited in the calorimeter cells, excluding those assigned to the outgoing positron. This 
reconstruction is sensitive to calorimeter noise effects. To improve it, especially at low masses, 
we removed isolated cells with energy below 140 (160) MeV for EMC (HAG) cells. The following 
formula was used for the mass reconstruction: 

4 



The coefficients Aq and Ai correct for the effects of energy loss in the inactive material and 
energy deposits below the threshold. Their values, Aq = 1.4 GeV and Ai = 1.21, were deter- 
mined from the Monte Carlo simulation so as to give the best estimate of the true invariant 
mass in diffractive events in the kinematic region of the analysis. In this region, the resolution 
in M^^'^ is approximately a{M^^'')/M^^'' ^ Q0%/^/M^. 

Diffractive events are characterized by a small value of Mx- However a simple selection based 
on M^"^ is not enough to select diffractive events because the non- diffractive background would 
be too high even at relatively low values of the mass. This is illustrated in Fig. ^3 where the 
In Mx distribution for the DIS sample with and without a forward rapidity gap requirement 
is presented along with a comparison with non- diffractive Monte Carlo expectations normalized 
to the integrated luminosity. Therefore, to reduce such background and extend the range of 
mass, a selection based on the presence of a forward rapidity gap is applied. The following cuts 
define the final sample for the analysis: 



Vmax — 1-8) where rjl^ax defined as the pseudorapidity of the most forward condensate 
with an energy above 400 MeV. A condensate is an energy deposit in contiguous cells 
above 200 MeV in the calorimeter. In Fig. |I]c the r/™^^ distribution for the full DIS data 
sample as well as for that with 7 < M^'^ < 25 GeV (see below) is presented along 
with a comparison with non-diffractive Monte Carlo expectations normalized to the total 
integrated luminosity. 

7 < M'^'^ < 25 GeV. The lower limit rejects events where the phase space for QCD 
radiation is limited. The upper cut is dictated mainly by the lack of statistics at high 
masses, a consequence of the rj^^ax requirement. 



W > 160 GeV. This cut helps to suppress the non-diffractive DIS background p2 



more than three condensates. This requirement ensures that an event plane for the 
hadronic final state X can be determined. 



A total of 2748 events are selected by these requirements. This sample is estimated to have 
a contamination from beam-gas background below 1%, based on the distribution of events in 
unpaired bunches. The background due to photoproduction is estimated to be less than 3% 
from Monte Carlo studies. Both backgrounds have been neglected in the analysis. 

The resulting reconstructed values of xjp are below 0.01, thus the selected sample is expected 
to be due to pomeron exchange 0, |3^ . The remaining non-diffractive DIS background varies 



from 5% at the lower end to 10% at the upper end of the mass range defined above (see 
Fig. P^). The background calculation has been checked by reweighting the exponential fall-off 
in the Monte Carlo sample in order to agree with estimates from fits to the data for In Mx ^'^'^ ^ 7 
with Mx'^ in GeV The two calculations agree within errors. This contribution will be 

subtracted on a bin-by-bin basis in the differential event shape distributions. 
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5 Monte Carlo simulation 



Monte Carlo (MC) event simulation is used to correct for detector acceptance and smearing 
effects. The detector simulation is based on the GEANT program |^6[ and incorporates our 
understanding of the detector and the trigger performance. 

Events from standard non-diffractive DIS processes with first order electroweak corrections 
were generated with the HERACLES 4.5 programme [0. It was interfaced to ARIADNE 



38| version 4.08 for modeling the QCD cascade according to the colour dipole model ||3^ that 



includes the boson gluon fusion diagram, denoted hereafter by CDMBGF. The fragmentation 
into hadrons was performed with the Lund fragmentation scheme as implemented in JETSET 

7.4 m. 



In order to model the diffractive hadronic final states we have considered various models for 



the pomeron. POMPYT is a Monte Carlo implementation of factorisable models for high 
energy diffractive processes, where, within the framework provided by PYTHIA ETI 
state proton emits a pomeron according to the Donnachie-Landshoff flux 



the initial 
E2|. The pomeron 

constituents take part in a hard scattering process with the virtual photon (see first diagram 
of Fig. 0a). The quark density in the pomeron is chosen to be of the type (3f{(3) oc (3{1 — (3). 
In POMPYT the final states come from the fragmentation of a qq pair, and QCD effects (e.g. 
QCD Compton, second diagram of Fig. |l]a) are approximated by the leading log approximation 
for parton showering. 



RAPGAP [0 represents an improvement over POMPYT in that a gluon density in the 
pomeron is also considered, thus giving rise to final states produced by a photon gluon mech- 
anism 7* (7 — > qq (see third diagram in Fig. |l]a). In the version we have used, the Streng 
parameterization for the pomeron flux [^] was used and the quark and gluon densities in the 
pomeron satisfy the momentum sum rule. Also, QCD effects are better modeled within the 

m 



colour dipole model 
shower approximation. 



5B[ than by the leading log parton 



as implemented in ARIADNE 
The colour dipole approximation turns out to be exact in the limit 
of massless quarks for the process of gluon emission from a qq pair. First order electroweak 
corrections are taken into account with the programme HERACLES 



A different approach has recently been proposed by VBLY [|1^]. Exact 0{as) calculations of 
two- and three-jet production in diffractive DIS are performed. The calculations assume that 
the pomeron is emitted at the proton vertex according to the Donnachie-Landshoff flux factor. 
The pomeron, assumed to be scalar, is considered to have a pointlike coupling to a quark pair, 
which may radiate a gluon, giving rise to qq or qqg final states, or to a gluon pair, leading to 
qqg final states through the fusion of one of the gluons with the photon at the lepton vertex 
(Fig. |T|a). Thus, the two free parameters in the model are the pomeron coupling to quark [gipqq) 
or gluon {gipgg) pairs. Two- and three-jet final states are generated using a Umin criterion]^ with 
the value set to 0.03. The value of the QCD scale parameter is chosen so as to reproduce 
the values for first order as as measured at PETRA. The Lund fragmentation scheme [P| 



as 



implemented with its default parameters in JETSET |^ is used. 



The RAPGAP and VBLY models are used for correction purposes, since POMPYT gives a 
much poorer description of the data, as discussed in the next section. 

^For any pair ij of partons, raf- > yminM"^ is required to avoid soft and coUinear divergences. 
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6 Results 



6.1 Global shape analysis 

In order to study the dependence of the multihadronic final states on Mx, we calculate in the 

n 

7*— iP rest frameQ the momentum tensor M^p = ''^^PiaPip, where the sum runs over the n 

i=l 

particles (condensates) in a final state and a, /3 = X, F, Z. Diagonalizing the symmetric tensor 
Majs yields three axes {k = 1, 2, 3) which give the orientation of the system together with 

n n 

three eigenvalues = ^^(Pi ■ n.k)'^, which can be normalized, Qk = Xk/ IPi^, such that 

1=1 1=1 
they satisfy the relation Qi + Q2 + Q3 = 1- By ordering the unit vectors Uk with increasing 

eigenvalues Qi < Q2 ^ Q3, one can define the variables sphericity S, the transverse momenta 
in (pt^^) and out (pTout) event plane as follows: 

n 

3 3 

o(<5i + Q2) = T^min^^^ {n = n3), 

" Erf 

i=l 



n 

i=l 



i=l 
1 " 



n — ' 

i=l 

where 723 is the unit vector along the principal axis, the so called sphericity axis, which minimizes 
the sum of the squared transverse momenta. The event plane is defined by n2 and na, while 
n\ defines the direction perpendicular to the event plane. Isotropic events are characterized 
by 5 ~ 1 and coUimated two-jets by 5 ~ 0. In models with a constant limited transverse 
momentum, the mean sphericity values are inversely proportional to the cm. energy. 



Alternatively, one can define the thrust |^ axis, denoted by the unit vector n-p, as the 
direction in space which maximizes the longitudinal momenta, with thrust then being defined 
as 

n 

' \Pi ■ n\ 



=1 

n 



T = max [n = nT^\ 



i=l 



Since thrust is a quantity linear in the momenta, it is less sensitive to QCD divergences. 
Isotropic events are characterized by T ~ 0.5 while two-jet events tend towards T ~ 1. 

Since the outgoing proton is not measured, the pomeron direction cannot be unambiguously determined. 
However, its transverse momentum is expected to be small. We therefore assume the emitted pomeron to be 
collinear with the incident proton, carrying a fraction xp oi its momentum. 
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In Figs. and 0, several observed differential distributions (i.e. uncorrected for detector and 
acceptance effects) are presented. Monte Carlo studies show that the detector acceptance and 
resolution do not significantly bias the shape of these distributions. 

The observed polar distribution ^|^'^^^ of the sphericity axis with respect to the virtual 
photon direction, for four different Mx'^ intervals, is shown in Fig. In e+e~ annihilation this 
distribution is found to be consistent with the form (1 + cos^Os) which is generally considered 
to give experimental evidence of the fermionic nature of quarks |^ . In contrast, the polar 



distribution of the sphericity axis measured in the cm. 7* — iP system is much more peaked, 
and, within the errors, little dependent on M^'^. This may be interpreted as a consequence 
of the t-channel quark propagator in Fig. |l]a for which the Born cross section is to a good 
approximation proportional to (1 — cos'^6s)~^- 

The observed sphericity and thrust distributions, in the four M^"^ intervals, are shown in 
Fig. 1^. As the cm. energy of the 7*— iP collision increases, the sphericity (thrust) distribution 
shrinks towards smaller (larger) values. This is a sign of coUimation, i.e. jet formation along 
the sphericity or thrust axes. 

The corrected mean sphericity (S), one minus mean thrust (1 — T) and the mean squared 
transverse momentum of the final state hadrons w.r.t. the sphericity axis (px), as well as its 
components in {pxij) and out {pr^^t) of the event plane, are shown in Fig. § as a function 
of Mx in the kinematic region 5 < < 135 GeV^ im < W < 250 GeV and r]^ax < 1-8. 
Here rjmax is defined as the pseudorapidity of the most forward generated particle with an 
energy above 400 MeV. These results are corrected for all detector effects, using Monte Carlo 
data passed through the detector simulation and data reconstruction program. The correction 
factor applied to the measured mean value, {X)'^^^^ , is given by the ratio (X)^^^ / (X)^^, 
where (X)^^^ and (X)^'^^ are the mean values from the generated and reconstructed Monte 
Carlo data respectively. (X)^^^ is evaluated at Mx and (X)^^^ at Mx'^. The correction factors 
using the RAP GAP and VBLY event generators deviate from unity typically by less than 5%. 
The errors quoted in this figure include, added in quadrature, systematic uncertainties due to i) 
model dependences in the correction procedure, ii) the determination of the boost to the 7*— iP 
cm. system using double angle or scattered lepton variables, iii) thresholds for the definition 
of condensates, iv) variations in the mass reconstruction procedure [Aq = and Ai = 1.5 in 
previous formula for Mx'^) and v) variation of the rj^^ax cut between 1.6 and 2.0. The systematic 
error dominates over the statistical. 

Correction to the whole r]max range is precluded by the fact that the correction becomes model 
dependent. In fact while the correction factor within the VBLY model is rjmax independent, in 
the RAPGAP model it drops to 0.6 at the highest mass bin. We interpret this as a consequence 
of the different behaviour of the pomeron remnant in the two models. In RAPGAP, the pomeron 
remnant, a quark or a gluon, follows the proton direction of flight, while in models based on 
pointlike couplings like VBLY, the partons at the pomeron vertex are produced at larger angles. 

The mean values (S) and (1 — T) measured in DIS events with a large rapidity gap at HERA 
(Figs. ^b) decrease with increasing Mx- This decrease is much slower than what would be 
expected in a model with constant hmited pt and opposite to isotropic phase space predictions 
^ . The mean transverse momentum with respect to the sphericity axis (Fig. ^) increases with 



increasing Mx- This growth is in fact mainly due to the transverse momentum component in 
the event plane, while that out of the event plane remains fairly constant (Fig. ^). Therefore, 
the events become planar. The observed features are similar to those in e~^e~ annihilation as 
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measured by the PLUTO ]T5[ and TASSO [||] collaborations at DORIS and PETRA (see Fig. |) 
at centre of mass energies comparable to the present range of Mx- Overall, the broadening 
effects measured in DIS final states with a large rapidity gap are stronger than those exhibited 
by e^e~ annihilation data. However, we note that the rjmax requirement precludes a model 
independent conclusion valid for the complete diffractive sample. 



6.2 Comparisons with models 

In the previous section we have presented evidence that the DIS large rapidity gap multihadronic 
final states become planar at high masses. POMPYT, for which broadening effects are given by 
the leading log approximation for gluon bremsstrahlung from the quark lines, does not describe 
the data. Not only does it fail to reproduce the distribution (not shown), but more 

important, the expected mean values for the event shape variables studied are much lower than 
the data (see Fig. ||). 



In the VBLY approach the two free parameters, gpqq and gpgg-, are directly related to 
the fraction of final state events coming from the gluon coupling as well as to the total diffractive 
cross section. This fraction was determined by fitting to a set of observed distributions such as 
sphericity, thrust, {pTout) ^^"^ |cos^5| as well as the (3 distributions shown in Fig. ||b, 

all of them plotted in the four mass intervals considered before. The parameters in RAP GAP 
have been tuned to describe the HI diffractive structure function F2 in and no additional 



tuning has been attempted in this analysis. 

The results of the VBLY fits (solid lines) and the expectations from RAPGAP (dashed lines) 
are shown in Figs. ^ 0, and ^. The global normalization in Figs. ^ and ^ is determined by 
the total number of events in the whole mass range 7 < M'^j^'^ < 25 GeV. Both models give a 
reasonable description of the data. Some discrepancies are observed in the | cos 6^5' | distribution 
in VBLY. Although in the VBLY model the pomeron coupling to quarks and gluons is pointlike, 
the (3 distributions, presented in Fig. are also reasonably described. The observed shape at 
small (3 is generated through the coupling of the pomeron to gluon pairs. 

RAPGAP with only a quark density in the pomeron (dotted line in Fig. ^) fails to describe 
the data. Similarly, VBLY with only the pomeron to quark pair coupling cannot describe the 
shapes of the sphericity, thrust and /? distributions (dotted line in Fig. |^ and |]). Note that 
the POMPYT prediction shown in Fig. ^ is similar to that of VBLY with only the quark pair 
coupling (not shown in Fig. 

The implementation of the hard gluon bremsstrahlung process in the MC models is insuffi- 
cient to account for the broadening effects observed in the data. Therefore, the need to include 
either the pomeron to gluon pairs coupling in VBLY or the gluon density in RAPGAP is very 
clear. The exact fraction of gluon-induced events is model dependent and it is sensitive to 
the i/min and pt cuts imposed on the LO matrix elements used. However, in both models a 
significant fraction of gluon-induced events is required to describe the data. In the kinematic 
range under study, this fraction is 50% for VBLY and 30% for RAPGAP. 
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7 Conclusions 



Wc have studied global event shapes in large rapidity gap NC DIS events at HERA in the 
kinematic range 5 < < 185 GeV^ 160 < < 250 GeV and ^rnax < 1-8. These events 
are generally interpreted as due to pomeron exchange. We have investigated the dependence 
of the event shape variables sphericity, thrust and transverse momenta squared in and out of 
the event plane with Mx, the cm. energy of the 7* — iP collision. The polar distribution of 
the sphericity axis is very peaked and little dependent on Mx- We find that with increasing 
Mx the large rapidity gap events become more collimated and planar. Broadening efi'ects in 
the DIS final states with a large rapidity gap call for a mechanism in addition to hard gluon 
bremsstrahlung. This can be achieved in models where the pomeron has a partonic structure 
by including a gluon density in the pomeron, as in RAPGAP, or through a pointlike coupling 
of the pomeron to quark and gluon pairs, as in VBLY. A significant gluon component of the 
pomeron is necessary for the models to describe the data. 
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Figure 1: a) Feynman diagrams for processes in the RAPGAP and VBLY models, b) InM^^^'^ 
distributions for the complete DIS as well as large rapidity gap samples, along with MC pre- 
dictions generated within the CDMBGF model and normalized to the integrated luminosity of 
the data, c) rj^^^ distributions for the complete DIS sample as well as for the restricted mass 
range 7 < M^'^ < 25 GeV, along with MC predictions as in b). 
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Figure 2: Measured polar distribution of the sphericity axis with respect to the virtual photon 
direction in the 7*— iP rest frame (a), and measured (3 distribution (b). The results of the VBLY 
model are shown as solid lines, with the quark contribution from this model shown as dotted 
lines. RAPGAP expectations are represented by the dashed lines. The global normalization of 
the predictions is adjusted so as to match the measured number of events in the whole mass 
range 7 < ^ < 25 GeV. 



14 



• ZEUS 1994 

- - RAPGAP 



a 



VBLY qq and gg couplings 
VBLY only qq coupling 




0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 

Sphericity Sphericity 

!/3 




0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 

Sphericity Sphericity 



b) 




0.5 0.6 0.7 0.8 0.9 1 0.5 0.6 0.7 0.8 0.9 1 

Thrust Thrust 




0.5 0.6 0.7 0.8 0.9 1 

Thrust 



0.5 0.6 0.7 0.8 0.9 1 

Thrust 



Figure 3: Measured sphericity (a), and thrust distributions (b) for large rapidity gap DIS events 
in four M'x'^ intervals. The results of the VBLY model are shown as solid lines, with the quark 
contribution from this model shown as dotted hues. RAPGAP expectations are represented by 
the dashed lines. The global normalization of the predictions is adjusted so as to match the 
measured number of events in the whole mass range 7 < M^^'^ < 25 GeV. 
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Figure 4: The mean sphericity, one minus mean thrust, mean squared transverse momenta w.r.t. 
the sphericity axis and its components in and out of the event plane, plotted as a function of 
Mx- Open dots and squares represent PLUTO and TASSO measurements at DORIS and 
PETRA, respectively. Black dots represent the ZEUS data corrected for detector effects to the 
kinematic region 5 GcV^ < < 185 GeV^ 160 GeV < < 250 GeV and r]^ax < 1-8. Note 
that the ZEUS data and the predictions from the models are for rjmax < 1-8. 
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